Method for detecting of battery condition

ABSTRACT

PROBLEM TO BE SOLVED: To provide a method for detecting a charged state of a battery, for evaluating a deterioration thereof due to each of reaction processes that individually have rates of reactions as different from therebetween, and for performing a detection of the charged state of the battery. 
     SOLUTION: The method for detecting the charged state of the battery according to the present invention comprises the steps of: measuring a voltage V mes  of the battery, an electric current I mes  thereof and a temperature T mes  thereof, and then inputting such the measured values, as a step S 1;  judging whether or not an absolute value of the measured electric current as the I mes  is smaller than a threshold value of the electric current as an I thre , as a step S 2;  estimating an OCV 20 hr  by making use of an SOC n−1  and an SOH n−1 , that are the values after charging and/or discharging at the last time, with reference to a stable OCV estimated formula, as a step S 4;  calculating a difference between the measured value of the voltage as the V mes  and the OCV 20 hr , and then saving such the calculated value, as a step S 5;  renewing a relaxation function as an F n (t) with corresponding to an amount of time as t, as processes from a step S 6  through a step S 19;  calculating an SOH n  at the step S 17  with making use of the F n (t) to be renewed; and calculating an SOC n  at the step S 19  with making use thereof.

CROSS-REFERENCE TO RELATED APPLICATIONS

This application claims priority from a Japanese patent application serial No. 2008-103925, filed on Apr. 11, 2008, the entire content of which is incorporated herein by reference.

BACKGROUND OF THE INVENTION

1. Technical Field

The present invention relates to a method for detecting of battery condition and a detection device for detecting the charged state of the battery.

2. Description of the Related Art

In recent years, there are made use of a lot of electric devices for a running of a motor vehicle, and then there becomes to be enhanced a significance of an electric power supply to be mounted thereon. In the past of twenty years to thirty years or longer before, a requirement for such the electric power supply to be mounted thereon is limited to a function, such as for starting an engine therein, for switching on an air conditioner therein or on a lamp thereon, or the like therefor. On the contrary, there becomes to be progressed on “by-wire” in recent years, and then there becomes to be controlled by making use of an electricity for a component part of a safety system, which is to be represented by an electric powered brake (EPB). Moreover, according to a conservation of energy, an emission regulation of carbon dioxide, or the like, as measures for an improvement of a fuel consumption, there becomes to be required a function of an idling stop at a period of stopping at such as an intersection or the like for a short period of time and to be required an ensuring of an ability for restarting such the engine therein. Thus, there becomes to be required a various kind of and a wide variety of functions for such the electric power supply therefor and for a battery therein, and then as corresponding thereto, there is desired an improvement on an accuracy for detecting a charged state of such the battery therein.

As following such the situation as mentioned above, to detect with an accuracy as higher regarding a state of charge (SOC) of a battery in particular becomes to be a technology as significant, for becoming involved in an stable operation of an electric device, such as the EPB or the like, for performing the running of the motor vehicle as safely and as comfortably, and then for realizing a social life based on the motor vehicle with taking into a consideration of an environment.

Under a condition in general for a battery as becoming to be stabled, an open circuit voltage (OCV) thereof and an SOC thereof have a relationship as corresponding one to one therebetween (refer to a symbol 81 in FIG. 12). However, the battery after performing a charging and/or a discharging is receiving each of effects of a generation and/or an annihilation reaction of ions due to an electrochemical reaction at a surface of an electric pole plate therein, and of a movement of the ions due to a diffusion and/or a convection of a battery electrolyte therein. And then thereby it takes time for convergent to the OCV to become stabled (approximately twenty hours for example). Moreover, in a case where there is such a change with the lapse of time, it cannot help but become apart from the correspondence as one to one for between the OCV and the SOC. Here, FIG. 13 and FIG. 14 are the graphs for individually showing one example of the change of the OCV with the lapse of time in a case where the SOC of the battery and a temperature thereof are individually constant respectively. Further, FIG. 13 shows that it takes time for the OCV thereof (refer to a symbol 82 therein) to become stabled as a predetermined value even in a case where the SOC thereof is constant. Furthermore, according to FIG. 14, there is shown each of variations of the OCV of the battery (refer to a symbol 83, a 84 and a 85 therein) in a case where a state of health (SOH) thereof is different from therebetween, and then the same shows that such the values are not convergent to the OCV as equivalent to therebetween if each of the SOH thereof is different from therebetween, even in a case where there is performed an adjustment for such as the SOC thereof, the temperature thereof, or the like, to be a condition as equivalent to therebetween respectively, and also in a case where there is designed a condition for the charging and/or for the discharging immediately before to be as equivalent to therebetween.

Thus, for evaluating the SOC thereof with reference to the OCV thereof by making use of the method for detecting the charged state thereof, there is not reflected any effect form the SOH thereof, if simply there is estimated the OCV thereof and then there is made use of the same, or if there is made use of only a hysteresis of the charging and/or of the discharging immediately before. And then in a case where the SOC thereof is evaluated with reference to the OCV thereof without reflecting thereto a condition of a deterioration regarding a battery, there become to occur a problem, such as that an accuracy of detecting the charged state thereof becomes to be worsened, or the like.

Here, as one example of the conventional technology, there is heretofore known a patent document 1. According to such the document, as a method for detecting the OCV of a secondary battery and the SOC of the same, there is made use of a transient response, which is in response to the hysteresis of the charging and/or of the discharging, for performing a compensation of the OCV thereof. Moreover, according to such the document, there is disclosed that the transient response thereof becomes to be varied in response to a period of time for the charging and/or for the discharging, and that there is made reference to such as the resistance component therein, the polarized component therein in response to the internal reaction of the battery, the diffusion velocity of the battery electrolyte therein, or the like.

[Patent Document 1] Japanese Patent Application Publication No. 2005-106615

However, according to the method for detecting thereof as disclosed in the patent document 1, it is not able to detect a deterioration thereof which is in response to a rate of reactions as slower, because it is not able to measure a transient response for a period of time as longer, though it is able to detect a deterioration thereof in response to a rate of reactions as faster at an inner portion of a battery in a case where a charging and/or a discharging thereof is performed for a period of time as shorter and then thereafter there is measured a transient response for a period of time as shorter. Moreover, it is necessary to perform the charging and/or the discharging thereof for a period of time as longer for measuring the transient response for the period of time as longer. However, an SOC thereof cannot help but vary as longer the period of time for the charging and/or the discharging thereof becomes to be, and then thereby the transient response for a shorter period thereafter cannot help but vary either. Thus, according to the method for detecting thereof as disclosed in the patent document 1, it is not able to catch each of the deteriorations in response to each of the rate of the reactions as different from therebetween at the inner portion of the battery, and then thereby to detect a charged state of the battery, such as a remaining capacity thereof or the like.

Further, according to a series of reaction processes, that is so called an electrochemical reaction, at an inner portion of a battery, there are effected to such the reaction processes not only by a generation and/or an annihilation reaction of ions to be generated in a vicinity of an electric pole plate (a rate of reactions as faster) but also by a diffusion velocity of ions in a battery electrolyte thereof (a rate of reactions as slower). Furthermore, according to such the reaction system as mentioned above, each of such the reaction processes, that individually have each of the rates of the reactions as different from therebetween, become to have an excessive effect as an error factor on an accuracy for detecting a charged state thereof.

Here, the present invention is presented for solving such the above mentioned problems, and an object thereof is to provide a method for detecting a charged state of a battery, wherein there is performed a detection of the charged state thereof by evaluating a deterioration due to each of reaction processes that has a rate of reactions as different from therebetween.

SUMMARY OF THE INVENTION

It is an object of the present invention to at least partially solve the problems in the conventional technology.

A method for detecting a charged state of a battery regarding a first aspect according to the present invention is characterized in that the method for detecting a charged state of a battery comprising the steps of: pre-forming a relaxation function F(t) for calculating a variation of an open circuit voltage (OCV) from the OCV at a time of a stable state regarding the battery after an elapsed time of t since the battery stops charging and/or discharging, as a function of a quantity regarding a predetermined state of the battery; measuring a variation of the OCV regarding the battery; optimizing the relaxation function F(t) with making use of the variation of the OCV to be measured; estimating the quantity of the state thereof by making use of the relaxation function F(t) to be optimized; and detecting the charged state of the battery with reference to the quantity of the state thereof to be estimated.

Although the invention has been described with respect to specific embodiments for a complete and clear disclosure, the appended claims are not to be thus limited but are to be construed as embodying all modifications and alternative constructions that may occur to one skilled in the art that fairly fall within the basic teaching herein set forth.

BRIEF DESCRIPTION OF THE DRAWINGS

FIG. 1 is a diagram showing a whole configuration regarding a detection system for detecting a charged stated by making use of the method for detecting the charged state of a battery regarding an embodiment according to the present invention.

FIG. 2 is a flow diagram explaining a judgment of a threshold by making use of the method for detecting the charged state thereof according to the present embodiment.

FIG. 3 is a flow diagram explaining an acquisition method for acquiring a value of voltage measurement by making use of the method for detecting the charged state thereof according to the present embodiment.

FIG. 4 is a flow diagram explaining a selection method for selecting an OCV_(20 hr) by making use of the method for detecting the charged state thereof according to the present embodiment.

FIG. 5 is a flow diagram explaining a calculation method for calculating a ΔV by making use of the method for detecting the charged state thereof according to the present embodiment.

FIG. 6 is a flow diagram explaining a calculation method for calculating a relaxation function F by making use of the method for detecting the charged state thereof according to the present embodiment.

FIG. 7 is a flow diagram explaining a calculation method for calculating the relaxation function F by making use of the method for detecting the charged state thereof according to the present embodiment.

FIG. 8 is a flow diagram explaining a calculation method for calculating an SOH and an SOC by making use of the method for detecting the charged state thereof according to the present embodiment.

FIG. 9 is a graph showing one example of a relaxation function and of a relaxation function for every rate of reactions.

FIG. 10 is a graph showing one example of a relaxation function for every rate of reactions and of a ratio therebetween.

FIG. 11 is a graph showing one example of a stable OCV estimated formula.

FIG. 12 is a graph showing a relationship between an open circuit voltage to be stabled and an SOC thereof.

FIG. 13 is a graph showing a variation of an OCV in a case where the SOC thereof is constant.

FIG. 14 is a graph showing a variation of an OCV in a case where the SOHs thereof are different from therebetween.

DESCRIPTION OF THE REFERENCE SYMBOLS

01: STORAGE BATTERY

02: DETECTION DEVICE FOR DETECTING CHARGED STATE

03: MOUNTED DATA CONTROL DEVICE

20: TEMPERATURE MEASUREMENT UNIT

21: VOLTAGE MEASUREMENT UNIT

22: AMPEROMETRY UNIT

23: STORAGE AREA (RAM)

24: FIXED STORAGE AREA (ROM)

25: ARITHMETIC AND LOGIC UNIT

26: JUDGMENT RESULT OUTPUT UNIT

27: MOUNTED STATE INPUT UNIT

28: TIMER

50: TRUE VALUE OF ΔV

51, 52, 53, 54: RELAXATION FUNCTION FOR EVERY RATE OF REACTIONS

55: RELAXATION FUNCTION

61, 62, 63: RELAXATION FUNCTION FOR EVERY RATE OF REACTIONS

64, 65: RATIO BETWEEN RELAXATION FUNCTIONS FOR EVERY RATE OF REACTIONS

70: STABLE OCV ESTIMATED FORMULA

81: SOC

82, 83, 84, 85: OCV

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS

A method for detecting a charged state of a battery regarding a preferred embodiment according to the present invention will be described in detail below, with reference to the drawings. Here, each of component parts that has a similar function is designated by making use of the similar symbol for simplifying to show in each of the figures and to describe thereof.

According to the method for detecting the charged state of the battery regarding the present invention, there is evaluated an SOH as an index of the state of health regarding the battery in response to a rate of reaction regarding the battery, and then by making use of such the SOH, there is evaluated an SOC thereof. Moreover, such the SOC has a relationship with an OCV thereof as one to one as shown in FIG. 12, and then by evaluating the OCV thereof, it becomes able to evaluate the SOC thereof as well. However, the relationship as shown in FIG. 12 is stood up for the OCV thereof in a case where the state of the battery is stabled. On the contrary thereto, the OCV thereof after charging and/or discharging has a variation as shown in FIG. 13. And then thereby it becomes necessary to evaluate such the SOC thereof by making use of the OCV thereof at a time therefor to become sufficiently stabled after the charging and/or the discharging.

Moreover, the variation of the OCV thereof after the charging and/or the discharging becomes to be changed due to the SOH thereof as shown in FIG. 14, and then thereby the OCV thereof at the time therefor to be stabled becomes to be changed due to the SOH thereof as well. Further, there is effected on a transient variation of the battery after the charging and/or the discharging not only by a matter having a rate of a reaction as faster, such as a generation reaction and/or an annihilation reaction of ions therein or the like, but also by a matter having a rate of a reaction as slower, such as a movement of a battery electrolyte therein or the like.

Therefore, even in a case where an elapsed time is shorter which is after charging and/or discharging, for evaluating an effect from the rate of the reaction as slower as well, and then for evaluating the SOC thereof as correctly, according to the method for detecting the charged state of the battery regarding the present invention, there is designed to be evaluated a variation of the SOH thereof for every rate of the reactions, and then by making use of such the SOH thereof there is designed to be compensated the OCV thereof.

Next, one embodiment according to the method for detecting the charged state of the battery regarding the present invention will be described in detail below.

At first, there is assumed an OCV_(S) at a period for becoming the OCV to be stabled after elapsing an amount of time as sufficiently since the charging and/or the discharging, and then there is assumed to be expressed for the relationship as one to one between the SOC thereof as shown in FIG. 12 and the OCV_(S) by making use of the following formula.

SOC=FS(OCV _(S)(SOC, SOH, T))

OCV _(S)(SOC, SOH, T)=lim(V _(mes)(t)).   (Formula 2),

Here, the above expressed “lim” means that an elapsed time as the t is performed to become infinite for a measured value of a voltage (OCV) as the V_(mes)(t) regarding a battery after elapsing the amount of time as the t since the charging and/or the discharging. Moreover, the OCV_(S) in the above expressed formula means that the same becomes to be varied due to the SOC thereof, the SOH thereof and the temperature as the T regarding the battery. Further, the V_(mes)(t) becomes to be varied due to the SOC thereof, the SOH thereof and the T regarding the battery as well. However, according to (Formula 2), there is expressed with making use of only the amount of time as the t when such the value becomes to be measured.

Next, in a case of a flooded type lead acid battery there becomes to be sufficiently smaller as not higher than 10 mV for a change with the lapse of time regarding the V_(mes)(t) per one hour at a time after elapsing twenty hours since an end of the charging and/or the discharging. Moreover, an error becomes to be as not higher than 0.1% for a degree of magnitude regarding the OCV thereof (approximately 12.9 V). Therefore, it becomes able to assume the V_(mes)(t) after elapsing twenty hours since the stop of the charging and/or the discharging to be as an OCV_(20 hr) in the following formula, and then there is designed to be made use thereof for the OCV_(S).

OCV _(20 hr) =V _(mes)(t=20 hr),

OCV _(S)(SOC, SOH, T)≈OCV _(20 hr)   (Formula 3).

Furthermore, it may be available to design a value for an amount of the elapsed time since the stop of the charging and/or the discharging as not equivalent to the twenty hours as well, with depending on a type of the battery.

Here, in a case where a variation from the OCV to be stabled regarding the V_(mes)(t) after stopping the charging and/or the discharging, that is to say, a variation of the OCV thereof, is assumed to be as ΔV(t), it becomes able to express the same by making use of the following formula.

ΔV(t)=V _(mes)(t)−OCV _(20 hr)   (Formula 4).

Moreover, such the variation of the voltage ΔV(t) has been treated as including all sorts of transient variations, with making use of a word as an “electrolytic polarization” according to a definition in the conventional electrochemistry. However, such the ΔV(t) is the variation of the voltage thereof which is generated due to a relaxation process in a period of time to become closing to the OCV thereof to be stabled, and then thereby such the value is influenced by a factor for the variation of the voltage thereof to be provided as below. For example, as a factor for the variation of the voltage thereof, there are provided such as a state of an electric pole plate therein, an concentration of ions in a vicinity of the electric pole plate, a solid state reaction therebetween, a solid-liquid reaction therebetween, a movement of the ions therein according to a precipitation of a battery electrolyte, to a convection thereof, or to a diffusion thereof, or the like. Thus, it may be considered regarding such the ΔV(t) as it is generated with incorporating the relaxation processes therein that individually have the rates of the reactions as different from therebetween.

Further, in a case where such the ΔV(t) is assumed to be expressed by making use of a function as an F(t), which is comprised of a polynomial expression having terms of an m as a natural number in response to a difference of the rates of the reactions therein, there becomes to be the following formula.

ΔV(t)=F(t)=f ₁(t)+f ₂(t)+ . . . +f _(m)(t)=Σf _(i)(t)   (Formula 5).

Moreover, according to the above expressed F(t) (the relaxation function), there is assumed that each of the terms as an f_(i)(t) (a relaxation function for every rate of reactions) express an extent of contribution thereby for the variation of the voltage thereof regarding each of the relaxation processes that the battery inherently has. And then each thereof individually becomes to be dependent on the state of health as the SOH of the battery, the state of the charging (concentration of the ions) as the SOC thereof, and the temperature as the T thereof respectively.

Further, there becomes to be designed a reference data to be stored at an inside of a detection system for detecting the charged state thereof before becoming to be connected between the battery and such the detection system for detecting the charged state thereof, with assuming as an initial state thereof in response to the battery to be connected therebetween, as expressed by the following formula.

SOC ⁰ =SOC ^(ref(0)) , SOH _(i) ⁰ =SOH _(i) ^(ref(0)) , OCV _(20 hr) ⁰ =OCV _(20 hr) ^(ref(0)).

Still further, for the initial state in which the battery and the detection system for detecting the charged state thereof are connected therebetween, there is designed to be made use of individual reference data for such the initial value respectively, as the measurement therefor at a number of times as an n equals to zero.

Still further, in a case where there are assumed to be as an F_(n)(t) and an f_(i) ^(n)(t) individually for the F(t) and each of the terms as the f_(i)(t) in (Formula 5) respectively, that are for expressing the variation of the OCV as the ΔV(t) after finishing the charging and/or the discharging at the number of times as the n (the n is an integer as not less than one) after becoming the battery and the detection system for detecting the charged state thereof to be connected therebetween and then becoming the initial value thereof to be set, it becomes able to calculate such the f_(i) ^(n)(t) with making use of an SOC thereof and an SOH thereof (that are assumed to be as an SOC^(n) and an SOH_(i) ^(n) respectively), that individually correspond to the rate of the reaction at a number of times as the i, by making use of the following formula.

f _(i) ^(n)(t)=f _(i) ^(ref)(t)(SOC ^(n) /SOC ^(ref))(SOH _(i) ^(n) /SOH _(i) ^(ref))g(T)   (Formula 6).

Here, the f_(i) ^(ref)(t), the SOC_(ref), and the SOH_(i) ^(ref) are the values that are pre-set as at the initial state thereof (a state as unused for example) respectively. Still further, the g(T) is a function that expresses a temperature dependency thereof.

Still further, in a case where there are assumed to be constant as independent of an amount of time for both of the temperature as the T and the SOC thereof according to (Formula 6), it becomes able to calculate the SOH_(i) ^(n) by making use of the following formula.

SOH _(i) ^(n)=(f _(i) ^(n)(t)/f _(i) ^(ref)(t))SOH _(i) ^(ref)   (Formula 7).

Still further, according to (Formula 7), an SOH^(n) in a total thereof, that there are summed up each of the SOH_(i) ^(n) due to each of the transient responses that individually have the rates of the reactions as different from therebetween, is comprised of each of components, such as expressed by making use of the following formula.

SOH ^(n)=(_(SOH) ₁ ^(n) , SOH ₂ ^(n) , . . . , SOH _(m) ^(n)), the m is a natural number   (Formula 8).

And then in a case where each of coefficients therefor is assumed to be as an A to an M for example, (Formula 8) becomes to be expressed by making use of the following formula.

$\begin{matrix} {\begin{matrix} {{SOH}^{n} = {{A\; {SOH}_{1}^{n}} + {BSOH}_{2}^{n} + \ldots + {M\; {SOH}_{m}^{n}}}} \\ {= {{{A\left( \frac{f_{1}^{n}(t)}{f_{1}^{ref}(t)} \right)}{SOH}_{1}^{ref}} +}} \\ {{{{B\left( {{f_{2}^{n}(t)}{f_{2}^{ref}(t)}} \right)}{SOH}_{2}^{ref}} + \ldots +}} \\ {{{{M\left( \frac{f_{m}^{n}(t)}{f_{m}^{ref}(t)} \right)}{SOH}_{m}^{ref}},}} \end{matrix}{{the}\mspace{20mu} m\mspace{14mu} {is}\mspace{14mu} {the}\mspace{14mu} {natural}\mspace{14mu} {{number}.}}} & \left( {{Formula}\mspace{14mu} 8^{\prime}} \right) \end{matrix}$

Still further, such the above expressed (Formula 8′) is just one example to be expressed in a form of the sum regarding the relationship among the SOH₁ to the SOH_(m) as shown in (Formula 8), and then regarding the total of the SOH_(i) ^(n), a form of the combination for each of the SOH₁ to the SOH_(m) is not limited to the form as expressed by (Formula 8′). Furthermore, it becomes able to detect a state of health regarding a battery with making use of such the SOH^(n).

However, in a case where an elapsed time is shorter since the end of the charging and/or the discharging at the number of times as the n, it is not able to evaluate the f_(i) ^(n)(t), that corresponds to the rate of the reaction as slower, and then thereby it becomes unable to renew the SOC^(n) and the SOH_(i) ^(n) therein. Therefore, before it becomes able to calculate the f_(i) ^(n)(t), the SOC^(n) and the SOH_(i) ^(n), that individually correspond to the rate of the reaction as slower, there are made use of the values at a time of the end of the charging and/or the discharging at the last time as an SOC^(n−1) and an SOH_(i) ^(n−1) in place of the SOC^(n) and of the SOH_(i) ^(n) respectively, and then there becomes to be made use of (Formula 6) with modifying into the following formula as approximately.

f _(i) ^(n)(t)=f _(i) ^(ref)(t)(SOC ^(n−1) /SOC _(i) ^(ref))(SOH _(i) ^(n−1) /SOH _(i) ^(ref))g(T)   (Formula 6′).

Moreover, according to the detection system for detecting the charged state thereof that it is able to apply (Formula 6′) thereto, regarding an measurement of a rate of relaxation regarding the F^(n)(t), it is able to perform the measurement thereof under a condition that the charging and/or the discharging of the battery is stopped. Moreover, in a case where there is always performed only an operation for the charging and/or the discharging thereof as not higher than a threshold therefor at all times thereof, it is possible to make use of the SOC^(n−1) therein for calculating the relaxation function as the F^(n)(t).

However, as assuming a condition of an operation with mounting on a motor vehicle, there becomes to be performed a charging and/or a discharging thereof according to the operation of the motor vehicle in a period of a measurement thereof at the number of times as the n after finishing the measurement thereof at the number of times as the n−1, and then there becomes to be required for compensating a ΔSOC (an added amount of the charging and/or the discharging) for the SOC^(n−1), as a changed portion of a charged amount in a period of the operation of such the motor vehicle. Further, in such the case thereof, there is assumed to be as the following formula:

SOC ^((n−1)′) =SOC ^(n−1) +ΔSOC.

f _(i) ^(n)(t)=f _(i) ^(ref)(t)(SOC ^((n−1)′) /SOC ^(ref))(SOH _(i) ^(n−1) /SOH _(i) ^(ref))g(T)   (Formula 6″).

Still further, there becomes to be renewed the f_(i) ^(n)(t) with making use of the SOH_(i) ^(n), that is calculated by making use of (Formula 7), and then such the renewed value becomes to be made use for performing the further calculation of the SOC_(i) ^(n) as expressed by the following formula.

f _(i) ^(n)(t)=f _(i) ^(ref)(t)(SOC _(i) ^(n−1) /SOC ^(ref))(SOH _(i) ^(n) /SOH _(i) ^(ref))g(T)   (Formula 6′″).

And then according to (Formula 4) and (Formula 6′″), it becomes able to calculate the OCV_(20 hr) by making use of the following formula.

OCV _(20 hr) =V _(mes)(t)−Σ[f _(i) ^(ref)(t)(SOC _(i) ^(n−1) /SOC ^(ref))(SOH _(i) ^(n) /SOH _(i) ^(ref))]g(T)   (Formula 9).

Furthermore, it becomes able to calculate the SOC^(n) by substituting such the OCV_(20 hr) into (Formula 2), and then it becomes able to make use thereof for detecting the charged state of the SOC thereof.

As described above, it becomes able to calculate the relaxation function as the f_(i) ^(n)(t) (i=1 to m, i and m are the natural numbers respectively) for every rate of the reactions as the number of the m pieces after finishing the charging and/or the discharging at the number of times as the n, with reference to each of the reference values as the f_(i) ^(ref)(t) (i=1 to m as the natural numbers respectively) as the number of the m pieces, which individually correspond to each of the rates of the reactions as the m types respectively, to the reference parameter for the state of health as the SOH_(i) ^(ref) (i=1 to m as the natural numbers respectively) as the number of the m pieces respectively, and to the reference parameter for the state of charge as the SOC_(i) ^(ref) (i=1 to m as the natural numbers respectively) as the number of the m pieces respectively. Thus, it becomes able to evaluate thereby the OCV thereof, the SOC thereof and the SOH thereof, that are individually reflected each of the states of health thereof respectively, that each thereof is in response to each of the rates of the reactions as different from therebetween respectively, and then thereby it becomes possible to perform a detection of the charged state thereof with an accuracy as higher.

Next, a method for detecting a charged state of a battery according to the present embodiment will be described in detail below, with making use of FIG. 1 through FIG. 7. Here, FIG. 1 through FIG. 7 are flow diagrams for explaining a flow of processes therefor by making use of the method for detecting the charged state of the battery according to the present embodiment.

Moreover, one example regarding a process of the method for detecting the charged state of a battery (01) according to the present invention will be described in detail below for a case of a battery which is mounted in a motor vehicle. Further, regarding a general view of the system therefor as shown in FIG. 1, a detection device for detecting the charged state (02) thereof is designed to comprise: a temperature measurement unit (20) for measuring the temperature of the battery (01); a voltage measurement unit (21) for measuring the voltage of the battery (01); an amperometry unit (22) for measuring the electric current of the battery (01); a storage region (RAM) (23) for recording temporarily each of measured values, which is measured by making use of each of the measurement unit (20 to 22) respectively; a fixed storage region (ROM) (24), in which a reference data is stored beforehand; an arithmetic and logic unit (25) for performing a detection and a judgment for the charged state thereof with reference to the data which is stored in the ROM (24); a judgment result output unit (26) toward an outside therefrom; a mounted state input unit (27) for being able to input data from a mounted data control device (03); and a timer (28) for being able to count an amount of time therefor.

Further, according to the present invention, there becomes to be executed a calculation operation of a quantity of a state thereof at a time for the battery (01) to be judged as the charging and/or the discharging becomes to be stopped. Here, there is shown in FIG. 2 regarding one example of a method for judging a stop of the charging and/or the discharging in the battery (01). Such the example is for a case, such as that there becomes to be judged by making use of the mounted data control device (03) that a motor vehicle is parking or is stopping, that there becomes to be input therein a data of putting on or of putting off the battery (01) and/or the detection device (02) for detecting the charged state thereof, that a value of an electric current, which is measured by making use of the amperometry unit (22) which is installed in the detection device (02) for detecting the charged state thereof, becomes to be not higher than a threshold for judging, which is recorded in the fixed storage region (24), or the like. Still further, it may be available to judge with making use of only the mounted state input unit (27), or to judge only the threshold of the electric current therein, or it may be available to judge with combining such the units in a free manner as well (here, such the above mentioned processes are defined to be as a threshold judgment process).

Still further, regarding the detection device (02) for detecting the charged state thereof, an initial state at a time therefor to be connected to the battery (01) is defined to be that a number of times as an n equals to zero, and then the number of times for the judgment thereof is defined to be as the n times at a time therefor to be judged by making use of the threshold judgment unit as shown in FIG. 2 that the charging and/or the discharging becomes to be stopped. Still further, a count rate by making use of the timer therein is defined to be as t_count=0.

Furthermore, according to the method for detecting the charged state of the battery regarding the present embodiment, there are assumed that SOC⁰=SOC^(ref(0)), SOH_(i) ⁰=SOH_(i) ^(ref(0)), OCV_(20 hr) ⁰=OCV_(20 hr) ^(ref(0)) regarding the SOC formula of (Formula 2), and then there becomes to be designed to save each of the values of the SOC of the battery, the SOH thereof and the OCV_(20 hr) thereof in the fixed storage region (24) beforehand, that each thereof individually corresponds to the initial state of the battery respectively.

Next, according to FIG. 3, there becomes to be performed a confirmation of the count rate of the timer therein with making use of a predetermined timing for the confirmation thereof, and then for a predetermined timing for a measurement thereof, at a time for the count rate of the timer therein (t_count) to become larger than the value of such the timing for the measurement thereof, there becomes to be measured a value of the voltage as the V_(mes) regarding the battery (01) by making use of the voltage measurement unit (21), and then there becomes to be determined that V_mes(t)=(t, V)=(t_count, V_mes) regarding a relationship between the amount of time and the voltage thereof at such the time thereof. Moreover, there becomes to be saved (V_mes(t)−OCV_(20 hr) _(—) temp) as ΔV(t)_temp (refer to (Formula 4); a process for obtaining the V_mes(t)).

Next, according to FIG. 4, there is shown a flow of a method for calculating the OCV_(20 hr) _(—) temp. Here, there are made use of the calculated values as the SOC^(n−1), and the SOH^(n−1) at the last time for the SOC and the SOH at a time of starting the calculation of the state for the number of times as the n. Moreover, regarding the temperature of the battery at the present, there is designed to be made use of a measured value T=T_n, which is obtained by making use of the temperature measurement unit (20). Further, as defining the value of the OCV after twenty hours to be as the OCV_(20 hr) _(—) temp, which is predicted with making use of such the above mentioned conditions, and then it becomes able to determine an OCV_(20 hr) ^(ref)(h) at the present, which is for a relational expression:

H(SOC _(—) j, SOH _(—) k, T _(—)1)=OCV _(20 hr) ^(ref)(h), (h, j, k, l therein are natural numbers respectively),

that has a reference data beforehand of which the OCV_(20 hr) ^(ref)(h) becomes to correspond to a condition of the h types that a plurality of the values of the SOC, of the values of the SOH and of the values of the T are combined therebetween. Thus, such the above mentioned process is the one of the embodiment regarding the stable OCV estimated formula. Furthermore, there becomes to be made use of such the value as:

OCV _(20 hd) _(—) temp=OCV _(20 hr) ^(ref)(h),

for the number of times as the n (a process for selecting the OCV_(20 hr) _(—) temp).

Next, according to FIG. 5, there becomes to be performed an addition of a data:

ΔV(t)_temp=V_mes(t)—OCV _(20 hr) _(—) temp,

into the temporary storage region (23) for every measurement thereof, with making use of the V_mes(t) which is calculated by making use of the measured value of the voltage according to FIG. 3, and with making use of the OCV_(20 hr) _(—) temp which is evaluated by making use of the quantities of the state as the SOC thereof, the SOH thereof and the T thereof according to FIG. 4.

Next, there becomes to be performed a fitting calculation for the ΔV(t)_temp, which is calculated according to FIG. 5, to be expressed as a sum of the relaxation functions for every rate of the reactions as the f_(i) ^(n)(t) (i=1 to m as the natural numbers respectively), that are individually formed as a form of function which is predetermined and as not less than two pieces (here there is assumed therefor to be as the m pieces).

In another word, there becomes to be performed an optimization of such the relaxation function for every rate of the reactions as the f_(i)n(t) for the data to be able to be expressed by making use of the ΔV(t)_temp.

Moreover, regarding a method for fitting thereof, there may be considered a variety thereof for calculating by making use of a regression calculation therefor, such as a method of least squares or the like. However, according to the present method for calculating such the ΔV(t)_temp, it cannot help but become the same to be as ΔV(20 hr)_temp=0, and then thereby the error thereof cannot help but become to be larger in a case where simply the regression calculation is performed therefor by making use of a sum of exponential functions.

Therefore, it is desirable to perform the following steps of: subtracting a gradient of a tangential line in a vicinity of the ΔV(20 hr)_temp=0; introducing a function for standing up that the ΔV(20 hr)_temp is larger than zero at all times thereof; and then performing a fitting by making use of the sum of the exponential functions for such the finite difference thereof.

Hereinafter, there is assumed (Formula 5) to be comprised of the four terms as expressed below for simplifying such the formula.

$\begin{matrix} \begin{matrix} {{F(t)} = {{f_{fast}(t)} + {f_{slow}(t)}}} \\ {= {\left( {{f_{{fast}\; 1}(t)} + {f_{{fast}\; 2}(t)}} \right) + {\left( {{f_{{slow}\; 1}(t)} + {f_{{slow}\; 2}(t)}} \right).}}} \end{matrix} & \left( {{Formula}\mspace{14mu} 10} \right) \end{matrix}$

Moreover, by expressing (Formula 10) as one embodiment for example to be as:

A function 1 for a rate of relaxation as faster:

F_fast1(t)=A exp(−B t̂C)   (Formula 10-1),

A function 2 for a rate of relaxation as faster:

F_fast2(t)=D exp(−E t̂F)   (Formula 10-2),

A function 1 for a rate of relaxation as slower:

F_slow1(t)=G exp(−H t̂I)   (Formula 10-3),

A function 2 for a rate of relaxation as slower:

F_slow2(t)=−a/72000 t+b   (Formula 10-4),

it becomes easier to perform a formation of a function to be optimized for the ΔV(t)_temp (here, each of the symbols “̂” therein means an exponentiation operator respectively).

Furthermore, it may be available to make use of such as another function to be further complicated than such the above expressed function, another function to be rather simplified, or the like, with depending on a condition of such as an operation speed of a sensor therein, a memory size thereof, an accuracy to be required for such the sensor therein, or the like.

Next, with making use of FIG. 6 and FIG. 7, a method for evaluating the F^(n)(t) by making use of the ΔV(t)_temp will be described in detail below. Here, there becomes to be set each of the coefficient for each of the functions as shown from (Formula 10-1) through (Formula 10-4) to be individually predominant as a fitting function individually for four of reference amount of times (10 sec, 1000 sec, 36000 sec, and 72000 sec) in such a time interval respectively.

And then at the time thereof, it becomes able to determine such the reference amount of times (10 sec, 1000 sec, 36000 sec, and 72000 sec) that are above mentioned as an example to be as in response to a time constant which is based on the rate of the relaxation regarding the rate of the reaction at the inner portion of the battery.

Moreover, it becomes able to change an amount of time that becomes to be such the reference therefor, with reference to an accuracy thereof and to a timing therefor that are required for a sensor therein, in response not only to the rate of the relaxation at the inner portion of the battery but also to such as a condition of a running by making use of an actual motor vehicle, a condition of a stopping of the same, or the like.

Further, regarding a method for determining such the above mentioned reference amount of time therefor, it may be available to make use of a timer at an inside of a sensor, or it may be available to make use of an amount of time that it is able to be obtained by making use of the mounted data input unit as well, which is represented by such as a car navigation or the like as a mounted information thereon, or it may be available to make use thereof with combining therebetween as well.

Still further, with making use of FIG. 6, the method for evaluating the F^(n)(t) by making use of the ΔV(t)_temp will be described in detail below, in a case where the elapsed time of the timer count therein is shorter than 20 hours.

Here, in a case where there is judged the amount of time as the t to be as shorter than a first reference amount of time (which is assumed to be as 10 seconds here), there becomes to be calculated the F^(n)(t) by making use of the following formula, with reference to the data as the F^(n−1)(t) after finishing the charging and/or the discharging at the last time thereof.

F ^(n)(t)=f _(fast1) ^(n−1)(t)+f _(fast2) ^(n−1)(t)+f _(slow1) ^(n−1)(t)+f _(slow2) ^(n−1)(t)   (Formula 11).

Hereinafter, as similar thereto, in a case where there is judged the amount of time as the t to be as not shorter than the first reference amount of time but shorter than a second reference amount of time (which is assumed to be as 1000 seconds here), there becomes to be calculated the F^(n)(t) by making use of the following formula, with reference to the data as the F^(n−1)(t) after finishing the charging and/or the discharging at the last time thereof and to the data as the latest.

F ^(n)(t)=f _(fast1) ^(n)(t)+f _(fast2) ^(n−1)(t)+f _(slow1) ^(n−1)(t)+f _(slow2) ^(n−1)(t)   (Formula 12).

Still further, in a case where there is judged the amount of time as the t to be as not shorter than the second reference amount of time but shorter than a third reference amount of time (which is assumed to be as 36000 seconds here), there becomes to be calculated the F^(n)(t) by making use of the following formula, with reference to the data as the F^(n−1)(t) after finishing the charging and/or the discharging at further the last time thereof and to the data as further the latest.

F^(n)(t)=f _(fast1) ^(n)(t)+f _(fast2) ^(n)(t)+f _(slow1) ^(n−1)(t)+f _(slow2) ^(n−1)(t)   (Formula 13).

Still further, in a case where there is judged the amount of time as the t to be as not shorter than the third reference amount of time but shorter than a fourth reference amount of time (which is assumed to be as 72000 seconds here), there becomes to be calculated the F^(n)(t) by making use of the following formula, with reference to the data as the F^(n−1)(t) after finishing the charging and/or the discharging at the most last time thereof and to the data as the most latest.

F ^(n)(t)=f _(fast1) ^(n)(t)+f _(fast2) ^(n)(t)+f _(slow1) ^(n)(t)+f _(slow2) ^(n−1)(t)   (Formula 14).

Still further, there is substituted the t=20 hours into the F^(n)(t) which is evaluated by making use of (Formula 11) through (Formula 14), and then there is assumed as:

OCV _(20 hr) ^(n) =F ^(n)(20 hr).

And then thereafter it becomes able to obtain the f_(i) ^(n)(t) and the OCV_(20 hr) ^(n) by making use of (Formula 11) through (Formula 14) and the following (Formula 15).

Furthermore, with making use of FIG. 7, in a case where there is judged the amount of time as the t to be as not shorter than the fourth reference time (which is assumed to be as the 20 hours for example), there becomes to be calculated the F^(n)(t) by making use of the following (Formula 15), with reference to the OCV_(20 hr) ^(n) which is evaluated by making use of the V_mes(20 hr) as the most latest and to the data of the ΔV(t)_temp that are recorded thereinto before, with assuming as:

ΔV(t)_(—) n=ΔV(t)_temp+OCV _(20 hr) _(—) temp−OCV _(20 hr) ^(n)   (Formula 14-2),

for such the ΔV(t)_n as not to be evaluated by making use of the F^(n)(20 hr), but as a form to be replaced with making use of (Formula 14-2).

F ^(n)(t)=f _(fast1) ^(n)(t)+f _(fast2) ^(n)(t)+f _(slow1) ^(n)(t)+f _(slow2) ^(n)(t)   (Formula 15).

Next, with making use of FIG. 8, there becomes to be calculated as the SOH_(i) ^(n), with making use of the fi(t), which becomes to be calculated according to the flow as shown in FIG. 6 through FIG. 7 and then the same becomes to be stored into the temporal storage region (23), and the reference data f_(i) ^(ref)(t) and the SOH_(i) ^(ref) that are recorded at the inside of the fixed storage region (24) beforehand.

Next, there becomes to be calculated as:

SOC ^(n) =SOC ^(ref) _(—) n,

by substituting each of the values into the relational expression of:

H(OCV _(20 hr) ^(n) , SOH ^(n) , T _(—) n)=SOC ^(ref) _(—) n,

which is stored in the fixed storage region (24), with making use of the OCV=OCV_(20 hr) ^(n), that is obtained according to FIG. 6 and FIG. 7, of the SOH=SOH^(n), that is obtained according to FIG. 8, and of the T=T_n, that is obtained according to FIG. 4.

Thus, it becomes able to perform the detection of the individual states regarding the state of charge regarding the battery and the state of health thereof, by comparing the SOC^(n) and the SOH^(n) to be calculated according to the above mentioned processes individually to each of the thresholds to be predetermined therefor respectively.

Next, regarding a comparison and judgment process, there becomes to be compared the SOC^(n) and the SOH^(n) individually to each of the thresholds to be predetermined which is pre-set therefor respectively, and then in a case where the SOC^(n) thereof is not larger than such the threshold therefor, or in a case where the SOH^(n) thereof is not smaller than such the threshold therefor, there becomes to be judged that the battery becomes to be at a deteriorated state. Furthermore, there may be a case due to a system design therefor where there becomes to be an inverse relationship between the large and small regarding the SOH thereof and the state of health regarding the battery, and then it is needless to say that there becomes to be judged that such the battery becomes to be at the deteriorated state in a case where the SOH^(n) thereof is not larger than such the threshold to be predetermined therefor.

Next, as one embodiment of the optimization of the relaxation function as the F(t) according to the method for detecting the charged state of the battery regarding the present embodiment, one example of a variation of each term in the relaxation function as the F(t) to be expressed by making use of (Formula 10) will be shown in FIG. 9.

FIG. 9 is a graph for showing the variation of the ΔV(t)(=F(t)) in a case where the horizontal axis therein is defined to be as an elapsed time since the end of the charging and/or the discharging, wherein each of the symbols from the 51 to the 54 designates each of the variation for each of the terms of (F_(fast1)(t), F_(fast2)(t), F_(slow1)(t), F_(slow2)(t)) in (Formula 10) respectively. Moreover, the symbol 50 designates the true value thereof, meanwhile, the symbol 55 designates the value for the F(t) which is calculated by making use of (Formula 10). Hence, with making use of such the above mentioned F(t) according to the present embodiment, there is shown in the figure that it becomes able to predict the ΔV(t) with an accuracy as higher.

Hereinafter, there is assumed to be as the SOH₁ in (Formula 7) for a degree of stratification regarding the rate of the reaction as slower (such as a diffusion of the battery electrolyte or the like), and then there is assumed to be calculated the SOH₁ ^(n) after finishing the charging and/or the discharging at the number of times as the n by making use of the following formula.

$\begin{matrix} \begin{matrix} {{SOH}_{1}^{n} = {{f_{slow}^{n}(t)}{f_{slow}^{ref}(t)}{SOH}_{1}^{ref}}} \\ {= \frac{\left( {{f_{{slow}\; 1}^{n}(t)} +_{{slow}\; 2}^{n}(t)} \right)}{\left( {{f_{{slow}\; 1}^{ref}(t)} + {f_{{slow}\; 2}^{ref}(t)}} \right){{SOH}_{1}^{ref}.}}} \end{matrix} & \left( {{Formula}\mspace{14mu} 16} \right) \end{matrix}$

Moreover, there is assumed in the above description that the f_(i) ^(n)(t) and the f_(i) ^(ref)(t) in (Formula 7) become to be calculated by making use of the sum of further two terms of (f_(slow1) ^(n)(t)+f_(slow2) ^(n)(t)) and (f_(slow1) ^(ref)(t)+f_(slow2) ^(ref)(t)).

As one example, with making use of the lead storage battery as the liquid type having the model number in size as 55D23 which is produced by FURUKAWA BATTERY, there are performed a cycle of the charging and the discharging from the unused state thereof as twenty times therefor, fifty times therefor and one hundred times therefor respectively, under a condition that an environmental temperature as 25° C. therefor and a depth of discharge (DOD) as 10% therefor. And then thereafter, with reference to the measured data after the cycle of the charging and the discharging as the twenty times therefor, after the cycle of the charging and the discharging as the fifty times therefor, and after the cycle of the charging and the discharging as the one hundred times therefor respectively, and with reference to the case of the cycle as the twenty times to be as the reference therefor, there becomes to be calculated the SOH₁ ^(n) by making use of the following formula.

SOH ₁ ⁵⁰=(f _(slow1) ⁵⁰(5 hr)+f _(slow2) ⁵⁰(5 hr))/(f _(slow1) ²⁰(5 hr)+f _(slow2) ²⁰(5 hr))SOH ₁ ²⁰   (Formula 17).

SOH ₁ ¹⁰⁰=(f _(slow1) ¹⁰⁰(5 hr)+f _(slow2) ¹⁰⁰(5 hr))/(f _(slow1) ²⁰(5 hr)+f _(slow2) ²⁰(5 hr))SOH ₁ ²⁰   (Formula 18).

Here, there is shown the value of (f_(slown)(t)/f_(slow) ²⁰(t)) in FIG. 10, that is calculated by making use of the measured data, wherein each of the symbols as the 61, the 62, the 63 individually designate the f_(slow) ²⁰(t), the f_(slow) ⁵⁰(t) and the f_(slow) ¹⁰⁰(t) respectively, and each of the symbols as the 64 and the 65 individually designate the (f_(slow) ⁵⁰(t)/f_(slow) ²⁰(t)) and the (f_(slow) ¹⁰⁰(t)/f_(slow) ²⁰(t)) respectively. Moreover, according to the same figure, at the point in time as the t=18000 seconds for example,

$\frac{\left( {{f_{{slow}\; 1}^{50}\left( {5\mspace{14mu} {hr}} \right)} + {f_{{slow}\; 2}^{50}\left( {5\mspace{14mu} {hr}} \right)}} \right)}{\left( {{f_{{slow}\; 1}^{20}\left( {5\mspace{14mu} {hr}} \right)} + {f_{{slow}\; 2}^{20}\left( {5\mspace{11mu} {hr}} \right)}} \right)} = {\frac{f_{slow}^{50}\left( {5\mspace{14mu} {hr}} \right)}{f_{slow}^{20}\left( {5\mspace{14mu} {hr}} \right)} = {1.52//}}$

As similar thereto, it becomes able to obtain:

$\frac{f_{slow}^{100}\left( {5\mspace{14mu} {hr}} \right)}{f_{slow}^{20}\left( {5\mspace{14mu} {hr}} \right)} = {1.63//}$

Thus, it becomes able to catch the variation of the state regarding the battery, with making use of the F(t) as the variation of the OCV, in response to the number of the cycles for the charging and the discharging.

Further, regarding the OCV after elapsing the time of 20 hours since the end of the charging and/or the discharging as one example, with corresponding to the number of the cycles for the charging and the discharging as 20, 50 and 100 respectively:

OCV _(20 hr) ²⁰=12.896 (V),

OCV _(20 hr) ⁵⁰=13.032 (V),

OCV _(20 hr) ¹⁰⁰=13.036 (V).

Still further, there is shown in FIG. 11 regarding a relationship between the above expressed (f_(slow) ^(n)(t)/f_(slow) ²⁰(t)) and the OCV_(20 hr) ²⁰. Furthermore, it becomes able to make use of the result as shown in FIG. 11 for a stable OCV estimated formula for estimating the OCV_(20 hr) regarding a battery of a type as similar thereto.

As mentioned in the above description, according to the present invention, it becomes able to provide the method for detecting the charged state of the battery, for performing the detection of the charged state of such the battery by evaluating the state of health thereof due to the reaction processes that individually have the rates of reactions as different from therebetween. Moreover, by performing the detection of the state of health as the SOH regarding the battery, it becomes possible to perform the detection of the state of charge as the SOC with the accuracy as higher. And then thereby it becomes possible to ensure the stable operation of the electric device, or it becomes possible to progress a prediction regarding any hazard therefor. Hence, there becomes to be obtained an advantage that it becomes able to maintain the running by motor vehicle in safe. Further, by improving the accuracy of the function for the idling stop, it becomes able to obtain another advantage as well, such as that it becomes able to reduce the load for the environment, or the like.

Furthermore, the description regarding each of the above mentioned embodiments is individually described for each example of the method for detecting the charged state of the battery according to the present invention respectively, and it is not limited thereto either. Regarding a detailed configuration, a detailed operation, or the like, according to the method for detecting the charged state of the battery regarding each of the above described embodiments, it is possible to modify properly within a scope that is not departing from the subject of the present

Thus, as described above, according to the present invention, there becomes able to provide the method for detecting the charged state of the battery, for performing the detection of the charged state of such the battery by evaluating the state of health thereof due to the reaction processes that individually have the rates of reactions as different from therebetween, and then thereby it becomes possible to perform the detection of the charged state thereof with the accuracy as higher. 

1. A method for detecting a charged state of a battery, comprising the steps of: pre-forming a relaxation function F(t) for calculating a variation of an open circuit voltage (OCV) of the battery after an elapsed time of t since the battery stops charging and/or discharging, as a function of a quantity regarding a predetermined state of the battery; measuring a variation of the OCV from the OCV at a time of a stable state regarding the battery; optimizing the relaxation function F(t) with making use of the variation of the OCV to be measured; estimating the quantity of the state thereof by making use of the relaxation function F(t) to be optimized; and detecting the charged state of the battery with reference to the quantity of the state thereof to be estimated.
 2. The method for detecting the charged state of the battery as defined in claim 1, wherein the relaxation function F(t) is pre-formed as a function of a temperature of the battery, and the temperature of the battery is measured and then the same is made used for the relaxation function F(t).
 3. The method for detecting the charged state of the battery as defined in claim 1, wherein a stable OCV estimated formula is pre-formed for calculating the OCV at the time of the stable state thereof, and the OCV at the time of the stable state thereof is calculated by making use of the stable OCV estimated formula and then a difference between the same and a measured value of the open circuit voltage of the battery is defined to be as the variation of the OCV.
 4. The method for detecting the charged state of the battery as defined in claim 3, wherein the OCV at the time of the stable state thereof is an OCV after an elapsed time of twenty hours since the battery stops charging and/or discharging.
 5. The method for detecting the charged state of the battery as defined in claim 1, wherein the quantity of the state thereof is a state of charge (SOC) and a state of health (SOH) regarding the battery.
 6. The method for detecting the charged state of the battery as defined in claim 5, wherein the battery becomes to be judged as abnormal, in a case where the SOC becomes to be as not higher than a first threshold to be pre-set, or in a case where the SOH becomes to be as not lower than a second threshold to be pre-set.
 7. The method for detecting the charged state of the battery as defined in claim 1, wherein the relaxation function F(t) is expressed as a sum of relaxation functions fi(t) (i=1 to m, i and m are the natural numbers respectively) for every rate of reactions as not less than two (assumed to be as the m terms), that is pre-formed with corresponding to each of the rates of the reactions regarding an inner portion of the battery, and each of the relaxation functions fi(t) (i=1 to m as the natural numbers respectively) for every rate of the reactions is optimized by dividing the measured value of the variation of the OCV into each of components for individually corresponding to each of the rate of the reactions.
 8. The method for detecting the charged state of the battery as defined in claim 7, comprising the additional steps of: measuring a voltage and an electric current of the battery; judging for the battery to be that becomes to be stopped charging and/or discharging, by making use of the electric current or of a predetermined signal regarding the stop charging and/or discharging; calculating a variation of the OCV from the measured value of the voltage, that corresponds to an elapsed time since the battery stops charging and/or discharging; optimizing each of the relaxation functions fi(t) for every rate of the reactions by making use of the variation of the OCV, that corresponds to each of the rates of the reactions which has a time constant as shorter than the elapsed time; making use of a relaxation function fi(t) as immediately before regarding the relaxation functions fi(t) for every rate of the reactions, that corresponds to each of the rates of the reactions which has a time constant as longer than the time constant; and estimating the quantity of the state by making use of the relaxation function fi(t) as immediately before and each of the relaxation functions fi(t) for every rate of the reactions to be optimized.
 9. The method for detecting the charged state of the battery as defined in claim 7, wherein the quantity of the state is calculated by estimating each of the quantity of the states for every rate of the reactions with making use of each of the relaxation functions fi(t) for every rate of the reactions, and by summing up the same.
 10. The method for detecting the charged state of the battery as defined in claim 9, wherein a relaxation function fin(t) for every rate of the reactions after finishing charging and/or discharging at a number of times as an n is expressed as fin(t)=firef(t)(SOCn/SOCref)(SOHin/SOHiref)g(T)   (Formula 1), in a case where the relaxation function fi(t) for every rate of the reactions, the SOC, and the SOH for every rate of reactions, at a reference condition for each thereof, are assumed to be as the firef(t), the SOCref, and the SOHiref respectively, and a dependency for a temperature thereof as the T is assumed to be as G(T) (here, the SOHin is the SOH for every rate of the reactions).
 11. The method for detecting the charged state of the battery as defined in claim 10, wherein the battery becomes to be judged that stops charging and/or discharging in a case where a value of the electric current is measured with which it becomes able to judge for a transient variation in the inner portion of the battery due to the charging and/or the discharging as substantially equivalent to a transient variation at a time of an open circuit therefor. 